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Abstract 



We present a general method for obtaining the quantum chromodynamical 
radiative corrections to the higher-twist (power-suppressed) contributions to 
inclusive deep-inelastic scattering in terms of hght-cone correlation functions 
of the fundamental fields of quantum chromodynamics. Using this procedure, 
we calculate the previously unknown 0{as) corrections to the twist-three 
part of the spin scaling function gxi^B, Q^)(= Qii^B, Q"^) + 92{xb, Q"^)) and 
the corresponding forward Compton amplitude STiv-iQ"^)- Expanding our 
result about the unphysical point xb = 00, we arrive at an operator product 
expansion of the nonlocal product of two electromagnetic current operators 
involving twist-two and -three operators valid to 0{as) for forward matrix 
elements. We find that the Wandzura-Wilczek relation between gi{xB-,Q'^) 
and the twist- two part of gxixB^Q'^) is respected in both the singlet and non- 
singlet sectors at this order, and argue its validity to all orders. The large- A'^c 
limit does not appreciably simplify the twist-three Wilson coefficients. 
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I. INTRODUCTION 



Deep-inelastic scattering (DIS) of leptons on the nucleon is a time-honored example of the 
success of perturbative quantum chromodynamics (PQCD) |l|]. The factorization formulae 
for the scaling functions Fi(xb,<5^) and QiixBiQ"^) (defined from the Bjorken limit of the 
structure functions Wiii'.Q'^) and Q^), respectively), augmented by the Dokshitzer- 

Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations for parton distributions [Q, 
describe the available DIS data collected over the last 30 years exceedingly well . Although 
the same formalism is believed to work for the so-called higher-twist structure functions 
1^, e.g. Wiiy^Q"^) and G2{y,Q'^), and their associated scaling functions, Fl^xb^Q"^) and 
g2{xB,Q'^), which contribute to physical observables down by powers of the hard momentum 
Q, there are few detailed studies of them in the literature beyond the tree level. [Actually, 
FiixsjQ"^) receives an as-suppressed twist-two contribution which we will ignore in this 
paper.] The QCD radiative corrections to Fl{xb,Q^) and g2{,XB,Q^) need be investigated 
as accurate data have recently been taken and more data will be available in the future 

i- 

In this paper, we present a general method for obtaining radiative corrections to the 
higher-twist parts of structure functions in inclusive DIS. This method is based on a gen- 
eralization of the tree-level formalism for higher-twist structure functions presented in Ref. 
1^, and has been outlined in a short paper of ours published earlier From this, we ex- 
tract a straightforward set of rules which can be used to calculate the amplitude of inclusive 
DIS at an arbitrary order in as and Aqcd/Q, where Aqcd is the QCD scale, in terms of 
light-cone correlation functions of quarks and gluons. As an illustration of our method, we 
present a detailed calculation of the radiative corrections to the lowest higher-twist observ- 
able, G2{i'iQ'^)- The partial results of this calculation have been communicated earlier in 
Refs. 1,0. 

In inclusive polarized deep-inelastic scattering off a nucleon target, one can measure the 
spin-dependent structure functions Gi^2{y, Q"^)- Gi{u, Q"^) is closely related to the spin struc- 
ture of the nucleon, and has been studied extensively, both experimentally and theoretically, 
in the last few decades |]TI|]. For quite a long period, the physics associated with G2{y, Q"^), 



in particular the scattering mechanism and its role in understanding the internal structure 
of the nucleon, has generated much debate in the literature |]T2|. In earlier investigations. 



attempts were made to interpret the structure function using the naive parton model and 



its various extensions |I3|. However, further studies showed that those results are incompat- 
ible with the factorization theorems of QCD [|1^]. Indeed, a QCD investigation shows that 
the scaling function g2{xB,Q'^) arises from the effects of parton transverse momentum and 
coherent parton scattering. The new information about the nucleon structure contained 
in G2{v-,Q'^) is found in its twist-three part [0, which reflects the quark-gluon correla- 



tions in the nucleon. In fact, this structure function represents the simplest manifestation 
of these correlations found in experiment. In addition, much of the information contained 
in G2{i',Q'^) has phenomenological interest. For example, the third moment of its scaling 
function, g2{xB,Q'^), is related to the response of the chromo-electric and -magnetic fields 
inside the nucleon to its spin polarization This information can help us to understand 



the role of the gluon fields in the nucleon, to test lattice QCD calculations, and to construct 
more accurate models of nucleon structure. 
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It is useful at this point to remind the reader of a fey key developments in studying 
the G2 structure function. Before QCD was invented, Burkhardt and Cottingham derived 
the super-convergence sum rule / du G2{i',Q^) = based on a dispersion relation and an 
assumption about asymptotic properties of the related Compton amplitude Jl^- In the 



context of QCD, a number of issues about G2{i', Q^) have been clarified in recent years. At 
tree level, g2{xB,Q^) is related to a seemingly simple quark distribution, qT{x,Q'^), within 
the nucleon [|TB|. However, as first pointed out by Shuryak and Vainshtein [^, its leading-log 
evolution reveals this simple form to be deceptive. In reality, one requires the introduction 
of a more general distribution to obtain a closed evolution equation [2U|. As found by 
Ali, Braun, and Hiller a somewhat surprising simplification occurs in the limit of large 
number of colors. In this limit, the leading order evolution of qrix, Q"^) becomes autonomous 
in the nonsinglet sector. However, it was later argued that this simplification does not occur 



at higher orders or in the singlet sector p2[. An analysis of the local operators associated 
with the tree-level expressions for the moments of g2{xB,Q'^) reveals a relation between its 
twist-two part and the scaling function gi{xBiQ^) |2^- This relation makes it possible to 
extract the twist-three part of g2{xB, Q"^) from the data. Experimental measurements of G2 



have been made by a number of collaborations ||2^. In addition, model calculations and 
lattice QCD yield interesting insights of the twist-three part of the structure function EH . 

This paper is organized as follows. In Section 2, a general discussion of the kinematics 
and structure functions associated with spin-dependent inclusive DIS is made. This section 
serves to familiarize the reader with some of our conventions and notations. Section 3 
is broken into five subsections dedicated to a systematic presentation of the method we 
use to obtain radiative corrections to the power-suppressed contributions to our amplitude. 
First, we summarize the well-known light-cone power counting techniques to the process 
at hand, which allows us to organize contributions to the amplitude with respect to their 
power suppression. Then, we discuss the role of longitudinal gluons in this process and the 
simplifications associated with the light-cone gauge. This is followed by a simple example 
in which the known tree-level result for 6*2(1^, Q^) is derived. It is here that the nonsinglet 
distributions in inclusive DIS with transverse polarization are introduced. Subsection D 
discusses the subtleties of renormalization and factorization relevant to radiative corrections. 
The last subsection presents a set of rules derived from this method which can be used to 
obtain the contribution to inclusive DIS at an arbitrary order ml/Q and ^^(Q^)- In Section 
4, we apply this method to the one-loop corrections of the results in Section 3C in the non- 
singlet and singlet sectors. These results were first presented in Refs. ^^e present 
communication, we go into somewhat more detail in explaining their derivation. An earlier 
analysis of the one- loop corrections in the singlet sector can be found in Ref . |^ . Although 
the preliminary amplitudes are the same, the final answer and its interpretation are quite 
different. A recent study of these corrections presents the same interpretation found 
here. There, one can also find discussions of the numerical significance of the results. Section 
5 presents the results of the last section in the form of an operator product expansion. In this 
form, it is easy to exploit the full space-time symmetries of our theory and isolate the new 
information provided by transverse scattering. This is one of the main results of this paper, 
as the twist-three Wilson coefficients have never been obtained at one-loop order before. In 
Section 6, we explain how to perform the twist separation without the use of the operator 
product expansion. Section 7 contains some concluding remarks and future possible work. 
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II. KINEMATICS 



To familiarize the reader with our conventions and notations, we review in this section 
a number of well-known facts about polarized deep- inelastic scattering (DIS). The process 
under consideration is inclusive lepton-nucleon scattering through one photon exchange. Al- 
though we have electron-proton scattering in mind, our results are applicable to processes 
involving any charged lepton and any hadron. Writing out the scattering amplitude, squar- 
ing, and summing over final states, one finds to leading order in Oem that the cross-section 
factors into two pieces : 

47rQ!eni t iivmit 

a oc -^L^^W^, ■ 

W^^{q,P,S) = ^Jd'ze'^-^{PS\[J^iz),rm\PS) , (1) 

where we have taken as the virtual photon momentum. J'^ is the electromagnetic current 
operator, 

nz)=Y.e,'^,{z)YM^) , (2) 
q 

and L^'^ is a tensor that depends only on the lepton spin and momentum. The ket \PS) 
represents a nucleon state of momentum P'^ and spin polarization S'^, which are constrained 
by = M^, S"^ = —1, and P ■ S = 0. The nontrivial kinematical invariants are = —q^ 
and u ^ q ■ P, which are positive in the physical domain. 

In polarized scattering, only the fiu antisymmetric part of W'^'^ contributes. Using dis- 
crete spacetime symmetries, one can express it in terms of two independent structure func- 
tions : 



M 



G,{u, Q^) + ±^{uS-q-SP)^ G,{u, 



(3) 



where [■ ■ ■] denotes antisymmetrization and e^'^°'^ is the Levi-Cevita tensor (gO^^a _ AH 
of the structural information on our nucleon state is contained in the dimensionless functions 
G'i,2('^, Q^)- Through the optical theorem, one can relate these structure functions to the 
invariant amplitudes contained in the forward Compton amplitude, 

T^"" = i j d^ze'''-' {PS\T {J^'{z)r{{))}\PS) , (4) 



as 



5i(z/, Q^) + ±-^[vS-q-SP)p S,{v, 



(5) 



G,{v,Q^) = ^\m[S,{v,Q^)] . (6) 

Owing to the time-ordered product in Eq. (H), the structure functions Si{h',Q'^) can in 
principle be calculated directly in Feynman-Dyson type perturbation theory. 

In the Bjorken limit, defined by taking Q^,z/ —>■ oo while the ratio xb = <5^/2z/ re- 
mains fixed, the physics involved in the structure functions becomes transparent because 
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the asymptotic freedom of QCD simplifies the scattering mechanism immensely. Well- 
defined relations emerge between the properties of the fundamental degrees of freedom — 
quarks and gluons — and the experimental data on the structure functions near this limit. 
The physical significance of the structure functions is made more apparant by introducing 
the scaling functions 



g,{xB,Q')^J^G,{u,Q') , 



(7) 
(8) 



which remain finite in the Bjorken limit. It is well-known that gi{xB, Q"^) contains informa- 
tion about the polarized quark and gluon distributions. From Eqs. ||), one finds that 
the structure function G2(z^, Q^) decouples from the scattering cross-section in the Bjorken 
limit. Hence its contribution is formally suppressed by a power of l/Q in comparison with 
Q^). At tree-level in QCD, g2{xB,Q^) is not related to a simple parton distribution. 
Rather, it contains information about quark-gluon correlations within the nucleon. The goal 
of this paper is to study the effect of QCD corrections on the relation between g2{xB,Q^) 
and these interesting correlations. 

At this point, it is convenient to choose the frame used most frequently in high-energy 
scattering. Taking the 3-momenta of and q'^ parallel and along the 2;-direction, we write 



M2 

where the basis vectors p'^ and are defined by 



p'^ = (l,0,0,l)A 



( is given by 



M2 



:i,0,0,-l)/2A 




(9) 
(10) 



(11) 
(12) 



(13) 



and A is an arbitrary scale reflecting our remaining freedom to boost along the 2;-axis. Note 



that this basis satisfies p 



n 



and p ■ n = 1. In the limit A ^ oo, one obtains 



the so-called infinite-momentum frame, in which Feynman's parton model was originally 
formulated p9i- the context of QCD, this limit can serve as the basis for light-cone power 



counting, which we will discuss in detail in the next section. 

Along with two independent transverse vectors, and n'^ form a complete basis for 
Lorentz vectors. Hence any vector, k^^, can be expressed in terms of them : 



A;^ = (n ■ k)p^ + {p ■ k)n^ + 
In particular, the spin polarization, S*^, has the form 



(14) 
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(15) 



Here, is an arbitrary vector in the 1- and 2-directions, normalized as Sj^ — —1. h 
and A represent the degree of longitudinal and transverse polarization, respectively, and are 
constrained by /i^ + = 1. In the limit A — > oo, the longitudinal polarization is leading 
since p'^ ~ A, whereas its transverse components are subleading since S'_l ~ 1. Taking 
A^ = 1, corresponding to complete transverse polarization, the hadron tensor simplifies to 



Prom the above, it is easy to see that one of the i^iy indices must be transverse while 
the other is longitudinal. More interestingly, it is grixBtQ'^) = gi{xB)Q'^) + g2{xB,Q'^) 
which naturally appears in DIS with a transversely polarized nucleon. For this reason, we 
will concern ourselves with this scaling function and the associated Compton amplitude 



Sriiy, Q^) = {u/M^)[Si(u, g2) + (u/M^)S2(iy, Q")], rather than g2(xB, Q') and S2{u, Q^). 



In this section, we present a general procedure to obtain perturbative QCD predictions 
for the structure functions in deep inelastic scattering, valid to all orders in Aqcd/Q and 
as{Q^)- For simplicity, we concern ourselves mainly with the forward Compton amplitudes. 
This procedure involves several steps. First, one needs to express the amplitudes in terms 
of general parton correlations functions and perturbative scattering amplitudes. Then a 
collinear expansion is used to simplify the classification of the contributions in terms of 
powers of Aqcd/Q- In the process, the parton correlation functions become light-cone 
correlations which have definite power counting properties. At a particular order in 1/Q, 
certain classes of processes differing in the number of longitudinal gluons can be resummed to 
form gauge links extrapolating between spacially separated physical fields in our correlations. 
This process of resummation is greatly simplified by the choice of light-cone gauge in our 
calculations. Once this choice is made, it becomes clear that only a finite number of gauge- 
invariant correlations can appear at each given order in 1/Q. The coefficients of these 
correlations are themselves power series in the strong coupling, as{Q'^)- As with most 
calculations in quantum field theory, higher order terms in these expansions contain infinities 
which must be understood and properly handled if one is to obtain a sensible result. After 
all of these steps are understood, one arrives at a simple set of rules which can be used to 
calculate the contributions to inclusive DIS at any order in 1/Q and as{Q^)- 

A. Collinear Expansion, Light-Cone Power Counting and Light-Cone Correlations 

As Q'^ —>■ oo, part of the deep- inelastic scattering process must be calculable in perturba- 
tion theory because of asymptotic freedom. On the other hand, the nucleon structure itself 
is surely non-perturbative. Therefore, the Compton amplitude T^'^ contains physics both at 
hard scales, 0{Q), and soft scales, O(Aqcd)- In the following discussion, we would like to 




gi{xB,Q'^) + g2{xB,Q'^) 



(16) 



III. THE FORMALISM 
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FIG. 1. Schematic diagram for Compton scattering in terms of quark and gluon scattering 
amplitudes M and their correlation functions F in the nucleon. 



separate the physics at these two scales. Although most of the discussion in this subsection 
may be found in the literature, the way it is presented here is new and general. 

As a first step of factorization, we consider the most general Feynman diagrams consisting 
of a hard (perturbative) part representing scattering of a set of quarks and gluons with 
the photon, and a soft (non-perturbative) part representing the Green's functions of these 
quarks and gluons in the presence of the nucleon, as shown in Fig. 1. Denote a collection 
of incoming (with respect to the hard scattering) quark momenta by {pi}, outgoing quark 
momenta by {p^} and gluon momenta by {ki}. The quark and gluon scattering amplitude 
can be expressed as 



= 1 j e^''<d'i{Q\T{j^{i)r{Q)---i,{p, 



>(p;)---A(fc,)---}|o) 



amp 



(17) 



where the subscript 'amp' indicates that the external quark and gluon legs are truncated 
and the Lorentz and Dirac indices on gluon and quark fields are left open. The quark and 
gluon Green's function in the presence of the nucleon is 



r(te}, {p;}, {h}) = {ps\t[-- i,{p,) ■ ■ ■ ij{p'^) ■ ■ ■ A{ke) ■ • ■}! P) . 
Therefore, the full Compton amplitudes can be written as 
T^''{q,p) = E / X{d'p^d%,d%M^^ (g, fe}, {p'^}, {h}) ® V (fe}, {p'^}, {h} 

all diagrams ijl 



(18) 



(19) 



where the coupling of Lorentz and color indices between M^*^ and F is implied. For simplicity, 
we work with bare (unrenormalized) fields and couplings until the end of the calculation. 

The contributions to T^^ can be classified in terms of powers oil/Q (loosely speaking, 
a twist-expansion). Unfortunately, however, 1/Q power counting is not simply dimensional 
counting. Conceptually, it is easy to understand because the transverse momentum of a 
parton, for instance, is less important than its longitudinal momentum in a hard scattering. 
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We must take these relevant facts into account systematically through a new power counting 
scheme. To facilitate this counting, we make a collinear expansion of the "hard amplitude" 
M'"'{q, {pi}, {p'j}, {kg}) by writing 

= Xipf" + Apf , 
p'/ = x'^p'' + Ap'^'' , (20) 
k'; = y,p^ + Ak^ , 

where Xi = Pi ■ n, and similarly for p'j and kg. Expanding all quark and gluon momenta 
about the collinear point Apf = Ap'j ^ = A/c^ = 0, we obtain 

M-to.{P.}.{p;}.{fc})= E ...J;,...,,,... 

[■ ■ ■ {Ap,r ■ ■ ■ i^PT'' ■ ■ ■ (^^^)*' ■■■][■■■ 5 - ■ -jM^^iq, {x.p}, {yep}) , (21) 

where we have omitted all the Lorentz indices which couple the partial derivatives with Api, 
Ap'-, and Akg. The leading term in the expansion can be interpreted as Compton scattering 
off of several collinear onshell partons with Feynman momentum fractions Xi = n ■ pi, etc. 
Higher order terms take offshellness and non-coUinear effects into account perturbatively. 

Rewriting all the quark and gluon momentum integrals in Eq. (|19[) in terms of these 
variables, we can integrate out Api, etc. Then the A-factors in Eq. (^if ) become spatial 
derivatives on the quark and gluon fields in T{{pi}, {p'j}, {ke}), which leads to 

T^''{q,p)= E E , — \ 

all diag.{m,},{n'^},{ti} ' ' ' ' ' ■ '^^ 

n / dx,dx',dye{- --d-- ■)M^''{q,p, {x,}, {x'j}, {yi})T^^t{{x^}, {x'^}, {ye}) . (22) 

ijk 

The new correlations Tmnt have field separations restricted to the light-cone ra^-direction in 
the coordinate space, and hence are called light-cone correlations: 

rmnt{{x.},{^-},{yi}) = n / ^^^e^(^-''-^-^-;+^^'^^) 

.^.^ J Zn Zn Ztt 

X {PS\ ■ ■ ■ d^'^^iXiu) ■ ■ ■ d"^^{njn) ■ ■ ■ d''A{uen) ■ ■ ■ \PS) , (23) 

where again we have omitted the Lorentz and color indices. 

Now, we are ready to introduce light-cone power counting. The hard part in Eq. (p2D 
contains only two momenta, and p'^, and the only dimensionful scalar in the absence 
of quark mass effects is Q^. [The effects of finite quark mass can be taken into account 
systematically, as explained in Subsection 3E.] On the other hand, the soft part depends 
on the scale Aqcd- Because there are Lorentz indices that couple the soft and hard parts, 
1/Q power counting is not determined entirely by the dimension counting of the soft and 
hard part alone. Instead we need a new type of power counting: light-cone power counting. 
Since the T^*^ has a fixed dimension, we can figure out 1/Q counting by studying behavior 
of either the hard or the soft part. Here we choose to examine the soft part. 
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The light-cone correlation functions in the collinear expansion involve QCD fields and 
their derivatives. Lorentz covariance requires that these objects be expressed in terms of p^, 
n^, = g^^ — p^n" —p^n^, e^^ = e^^^^PaUfj, and ^j^^. Different 1/Q powers result when 
these vectors are contracted with the hard partonic scattering amphtudes. Note that the 
contractions of Lorentz structures never produce any soft dimension because of the choice of 
the collinear coordinates. The balancing soft dimension in each term comes from the scalar 
coefficients of the expansion, taking into account [p] = 1, [n] = — 1 and [g] = [e] = [S±] = 0. 
Because the dimension of the Compton amplitude is fixed, the soft dimensions of the scalar 
coefficients determine the associated hard dimensions. Obviously, the terms with the lowest 
soft dimension dominate the contribution. 

For example, if the soft and hard parts are connected by two-quark lines, we have the 
following light-cone correlation function : 



/ 



^ e^'-{PS\UO)M^n)\PS) = 2 [j^^j^A^) + j'^^nMx) + • • ■] (24) 



By dimensional counting, [q{x)] = and [/4(a:)] = 2. Therefore, the q{x) contribution to 
the hard scattering is leading and the /^{x) contribution is suppressed by a factor of 

relative to q{x). 

The relation between counting the soft dimension and the use of the infinite-momentum 
frame is now easy to see. Since p^ ~ A, the contribution of the associated term is leading 
in the limit of A — >^ oo. On the other hand, n^^ ~ A~-^, so contributions associated with this 
vector are doubly suppressed. The Lorentz structures spanning transverse dimensions are 
0(A°) as A — oo. Hence, their contributions are suppresed relative to the p'^ terms, but 
enhanced related to the type. 

This fact allows us to separately identify a certain power suppression for each field and 
derivative in the correlation function. For gluon fields, one has 

~ (l/g)° , ~ (l/QY , A- ~ (l/Qf . (25) 

As required by the form of the covariant derivative, partial derivatives obey the same count- 
ing rules as gauge fields. Now, we see why our power counting arguments cannot be made 
until after the collinear expansion : subleading momentum components alter the form of the 
correlation functions, leading to an extra suppression. Until this behavior is made explicit, 
our correlation functions do not have a definite power counting. An analysis of the quark 
bilinears, taking into account the fact that fermion fields have mass dimension 3/2, leads to 

_V^7+(1,75)V' ~ (1/Q)' , ^7^(1,75)^' ~ (1/Q)' , ^7"(1,75)V' ~ (1/Q)' , 

i^ii, 75)V' ~ (1/Q)' , ~ (1/Q)' , V'^^^V ~ (1/Q)' , (26) 

where 7^ = (7°±7^)/\/2 and (1,75) means either 1 or 75. These relations can be simplified 
by breaking the quark field into two parts : 

^ = + , V'± = ^7VV' , (27) 



and writing simply 



^+ ~ (1/Q)i , V'- ~ (1/Q)' ■ (28) 
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The suppression associated with an arbitrary correlation function is simply the sum of the 
suppressions associated with its constituent fields and derivatives minus two for the external 
nucleon states. 

A general light-cone correlation F is the nucleon matrix element of an operator with 
(j) fields, for (p = A'^ , A^, A~ ,ip+,ilj-, and rig partial derivatives, for 9 = d±,d^. According 
to the above discussion, its soft dimension is given by 

noj^ + riA^ + + 2{nd~ + ua- + ) + ■ nA+ - 2 . (29) 

Armed with this formula, we can immediately determine the order at which a particular 
correlation function enters our amplitudes. For instance, the leading contribution to DIS 
comes from correlations with either two ip^ or two Aj_ correlations with arbitrary number 
of A+. 

For polarized DIS with transverse polarization, the leading contribution is of order 1/Q. 
The possible field combinations at this order are 4':f4'±, A_A±, A±A±A±, ip+A±ip+, ip+d±ip+ 
and A±d±A±, each appearing with an arbitrary number of A'^^s. For each of these correla- 
tions, we need to consider all possible Feynman diagrams with corresponding external legs. 
The momenta entering these legs are all coUinear. 

B. Longitudinal Gluons, Gauge Invariance and the Light-Cone Gauge 

As we have seen in the last subsection, the longitudinal gauge potential, A'^, in light-cone 
correlations does not lead to any suppression in its contribution to scattering amplitudes. 
At any given order in 1/Q, an infinite set of correlations with increasing number of A'^^s 
contribute. A related problem is that of gauge-invariance. The correlations that naturally 
appear in the collinear expansion are not manifestly gauge invariant. In particular, the 
spatially separated fields do not form gauge-invariant operators without appropriate gauge 
links. The appearance of gauge potentials also poses a gauge-invariance problem since 
these objects do no transform covariantly under the action of the gauge group. As we 
have mentioned before, the parton scattering amplitudes M'^" are calculated with external 
collinear and on-shell parton states after collinear expansion. As such, they are separately 
gauge invariant. Since the Compton amplitudes are gauge invariant order by order in 1/Q, 
all correlations at a given order must combine in such way to yield gauge invariant parton 
distributions. Gluon potentials and partial derivatives do not transform covariantly under 
the gauge group, so they must combine to form either covariant derivatives or field strength 
tensors. 

In this subsection, we argue that the infinite number of correlations must be combined 
into a finite number of gauge-invariant parton distributions. In particular, the correlations 
with increasing numbers of ^''''s form gauge- invariant structures with gauge links connecting 
fields at separate spacetime points. This conclusion is, of course, well-known in the case of 
leading twist We assert that it is in fact true order by order in 1/Q. 

Let us recall the argument of gauge-link formation in the leading twist case. Consider 
a class of diagrams with a number of A'^ gluons connecting the soft and hard parts. As 
we stated above, the hard part is already gauge invariant and can be calculated in any 
gauge one likes. Using Ward identities derived from the gauge symmetry, one factorizes 
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(b) 



(e) 



FIG. 2. A set of tree-level twist-two contributions to DIS, differing in the number of longitudinal 
gluons participating in hard scattering. 



all the longitudinal gluon lines connected to the hard scattering blob onto eikonal lines. 
These eikonal lines represent the gauge links. To prove this fact, one refers to the proof of 
factorization by Collins, Soper and Sterman in Ref. 0]. At leading order in l/Q and as{Q'^), 
this process is especially easy to see. Collecting the contributions from the diagrams in Fig. 
2, one obtains 



dxdX 
2n 



PS 



+ 



^+(0)7'^ 
dyidjjLi 



27r 



+ 



dyid/ii (i?/2(i/i 



2tt 



271 



2 g«A'lJ/l+«A'2l/2 



X + yi) ^+ ^ 
i 



-igA^{fiin)Y 



X p+ 



X 



{x + yi) i 



a; + 2/1 + V'l) ^+ i 
-igA'^ifiin)'^ 



X 



igA'^{lj,2n)Y 
7>+(An)|P5) 



(30) 



For illustrative purposes, we take fiu in the transverse dimensions and simplify the Dirac 
algebra. Giving xb = Q^f^p-Q a small negative imaginary part and performing the integrals 
over the yi, we obtain 

^2 



rplMU 



dxdX 



AXx 



X 



27r 





X 



Xb 



PS 



1 + / dpii 



A 



/■O rill 

-igA^{lJin)) + / dfii {-igA'^{fiin)) / dfi2 {-igA'^{fi2n)) + ■ ■ 
Jx J\ 



X 7XAn)|PS') + (/i ^ z^, g -g) 



(31) 



The quantity in brackets is precisely the gauge link required to make this distribution gauge- 
invariant. 

At higher orders in the 1/Q expansion, we have more fields and this process of resumma- 
tion becomes more complicated. However, at any given order in 1/Q, there is only a finite 
number of A± and A~ according to power counting. These potentials require only a finite 
number of A+'s and partial derivatives to form gauge covariant quantities. Therefore the 
infinite number of remaining A'^^s must form gauge- covariant objects by themselves. One 
can easily see by inspection that this can only be achieved if they form gauge links. Once 
done, we are left with a finite number of correlations with gauge links. 
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After an infinite number of longitudinal gluons have been resummed to form the gauge 
links, the remainder must combine with the partial derivatives and the A± and fields to 
form fully gauge-invariant distributions. This means that they must organize into covariant 
derivatives and field strength tensors. Although difficult to prove, this fact is a direct 
consequence of gauge invariance. It can be used as a check of any explicit calculation. 

In a practical calculation, the presence of the longitudinal gluons represents an unnces- 
sary complication. The standard approach to avoiding this is to make the choice of light-cone 
gauge ^4+ = 0. Although a gauge choice here is not required since we are not making per- 
turbative calculations of the correlations, it represents a tremendous simplification. In this 
gauge, no correlations involving ^4+ are ever needed and the gauge links simply reduce to 
unity. The only complication comes at the end of the calculation, when we would like to 
organize the correlations in the axial gauge into manifestly gauge invariant ones. 

As a last point in this subsection, we mention that not all the correlations that one 
can write down according to the light-cone power counting above are independent. The 
so-called 'bad' field components, ■0- and A~, are related to the 'good' components, ■0+ and 
A±, through the QCD equations of motion: 

flipi) =i^iD^ip+ flip^ip = 2iD+il)_+ flip±ip+ = , (32) 

where we have introduced the covariant derivative 



DaiiJ^n) = da + igAainn) . (33) 

In the light-cone gauge, we can solve this equation for the 'bad' quark field component in 
terms of the 'good' one: 

^_(An) = \f dC fli]/)^{CnyMCn) . (34) 

Similarly, one has 

A-{Xn) ^-fdC f dC [Df{Cn)d+Ai{Cn) + gi^+{C,'n) if^i^^iC'ri)] . (35) 

from the gluon equation of motion. The above equations are consistent with the light-cone 
power counting of the previous subsection. In both cases, we see that the 'bad' components, 
-0- and A~ , can be eliminated in favor of either two 'good' fields, ■0+ or ^4^, or one 'good' 
field and a transverse momentum insertion. This follows from the fact that in light-cone 
quantization, 'bad' field components do not propagate freely and hence cannot be directly 
interpreted in Feynman-Dyson type perturbation theory. As a consequence, we can choose 
independent correlations without inclusion of any "bad" fields. 



C. Example: St in the Leading Order 

As an example of formalism presented above, we consider the leading contribution to 
S'r(z/, Q^). This example is useful not only as a simple illustration of the procedure, but also 
to familiarize the reader with some of the distributions that appear in DIS with transverse 
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(a) 



(b) 



(c) 



+ 



FIG. 3. Tree-level Feynman diagrams for St- The pluses (+) and minuses (-) refer to the 
light-cone projections of the Dirac fields. 



polarization. To simplify the discussion, we contract one of the photon polarizations with 
n^. This leads to 



It 



MS, 



/3 



(36) 



where fi is transverse. 

In the light-cone gauge, there are two independent correlations at leading order in as{Q'^). 
The first involves the good components of two quark fields and a transverse partial derivative: 



qqB 



27r27r* 



{PS\MO)tdlM>^n)\PS) 



5{x-y) I ^e^'^^{PS\^P+{0)^^lM>^n)\PS) 
Zn 



(37) 



where the Dirac and color indices on quark fields are open and a is transverse. The second 
correlation contains two quark and one gluon field : 

^^e^'^e^^(^-^)(P5|^+(0)(-^7^)Al(/in)^+(An)|P5) . (38) 

The leading-order result for 5't(z^, Q"^) must be expressible in terms of these correlations 
alone. 

The tree level Feynman diagrams for Af^"*" are shown in Fig. 3. The total contributions 
to T^^ can be expressed as 



T^^-^ = ^ fdxdyf'-^^ 



dX 



1 



y-XB 



■7^ fii]p^{^n) 



1 



X + Xb 



X l^PS 

/ dx l-e^'^(PS 
2uJ J 27r 



PS 
PS 



PS 



(39) 

(40) 
(41) 
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The special form of the kernels in contribution (|39|) allows us to "simplify" this contribution 
using the equation of motion. It then combines with the others to give the full result 



— a/ 



Xb) tt- dx 



2u 



1 



X ~ Xb X + Xb 



X / — e 

2n 



^(0)7a75^(Ari) 



PS 



(42) 



As a consequence of parity and time-reversal invariance, one can isolate a scalar distribution, 

?(0)7a75^(An) PS) = {MS)o^qTB{x) , (43) 



dX . 
— e 
in 



iXx /pg 



and write 



dx 



X — Xb X + Xb 



qTB[X) 



(44) 



at leading order in the strong coupling. Recall that this distribution, as it appears here, is 
not renormalized. However, since the difference appears only at higher orders in ^^(Q^), we 
can replace it with a renormalized distribution qxix, Q^). [Note that qrix, Q^) can be made 
manifestly gauge-invariant by re-introducing the fields in the form of a light-cone gauge 
link connecting the quark fields.] 

It is now a simple matter to determine the scaling function Qri^BiQ'^) in terms of 
qrix^Q^) : 



qrixB, Q^) + qri-XB, Q' 



(45) 



This extremely simple result is quite misleading. It seems to imply that only one function 
of a single variable is required to predict the 0{1/Q) corrections to DIS. In analogy with 
the leading 0{1) result, one might assume that higher order corrections will do nothing 
but modify the form of the coefficient function. Furthermore, one might expect that the 
variation of this new distribution is fixed from the distribution itself. Unfortunately, 
neither of these suspicions is correct. The fact that our leading result can be expressed 
so simply in terms of qxiXyQ"^) is due exclusively to the form of the coefficient function in 
contribution (|39|). There is no reason to believe that this simple form will persist at higher 
orders. 

To reveal the full content of QrixBiQ'^) at this order, we use the equation of motion to 
eliminate the "bad" fields. The result can be expressed in terms of the single gauge-invariant 
correlation 



KqB^x, y) = T'^gBix, y) + r^ggB(2;, y) 



27r 27r 



{PS\ilj{Q)iD1{fm)^{Xn)\PS) 



(46) 



where color indices have been contracted but the Lorentz indices are open. The only invariant 
correlation KqB{x,y) appears in the contraction 
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J 4:71 4:71 \ 



PS)^ ie-+''mSaKgB{x,y) , (47) 



f ^ ^ gi^jz-o.) /pg ^(Q)^M ^ip^(ptn)tlj{Xn) PS) = -ie-^'^'MSaKgeiy^x) , (48) 



47r 47r 

of these latter indices. Since 



r ^rf^^axgiM(2/-a.) /ps ^(0)(-i^)F+°(/xn)7„ ^7'^V'(An) 
i 47r 47r ^ 



P5 



= ie-+"'^M5„(|/-a;)ir,B(a;,|/) , 



(49) 



every gauge-invariant way to include two 'good' quark fields and either a 'good' gluon field 
or a transverse partial derivative can be accomodated by KgB{x,y). Furthermore, the first 
moment of KgB{x,y) with respect to y generates the 'bad' component of the quark field 
through the equation of motion. Hence all of the distributions relevant to 0{1/Q) DIS are 
contained in KqB{x,y). This fact allows us to assert that all 0{1/Q) corrections to DIS in 
the nonsinglet sector can be expressed in terms of KqB{x,y). In the singlet sector, one must 
introduce one more distribution to describe pure gluonic effects. This is done in Section 4. 
In terms of KqB{x, y), our leading correction has the form 



n J X \X — 



Xb X + Xb 



Kg{x,y,Q'') 



(50) 



where we have once again ignored the difference between renormalizcd and unrenormalized 
distributions. In the remainder of this paper, we will be concerned with determining the 
0{as) corrections to this result. 



D. Factorization of Infrared Singularities and Radiative Corrections 

Up to this point, we haven't said much about the perturbative coefficients M^" of the 
correlation functions other than that they are calculated in terms of collinear and on-shell 
external parton states and that they are gauge invariant. They are, of course, perturbative 
expansions in the strong coupling, as{Q^)- Hence, the final result for T^^ is a double 
expansion. As Q becomes large, both expansion parameters vanish. Loosely speaking, the 
1/Q expansion is associated with the number of fields in the soft part of the process, while 
the coupling expansion is associated with those in the hard part. In this subsection, we 
would like to elucidate the nature of the latter expansion. 

Because the partons are on-shell, Mf^ can be viewed as a parton scattering S-matrix 
element. In principle, one must multiply by parton wave function renormalization factors 
to get the proper S-matrix element. However, working in d = 4 — e dimensions, the absence 
of a physical scale at the massless poles of the quark and gluon propagators reduces these 
contributions to unity. Subleading terms in the expansion of Mf*^ are parton scattering 
S-matrix elements with insertions of certain vertices associated with powers of kj- and k~{= 
p-kj). For instance, with one power of quark momentum kj-, the subleading term is calculated 
with one insertion of the vector vertex to one of the quark propagators. Because S-matrix 
elements and their relatives are gauge invariant, one may choose any gauge for the internal 
gluon propagators in M-^'^ . We use Feynman's choice in our calculation. 
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As we have argued above, M/"^ is ultraviolet finite because the on-shell wave function 
renormalization is trivial in dimensional regularization. Nevertheless, due to the massless 
on-shell external states, M/*^ has infrared divergences showing up as 1/e poles. To under- 
stand the origin of the infrared divergences, we consider a vertex correction to the tree-level 
Compton scattering diagram. Let the gluon momentum be k and one of the quark propaga- 
tors k + xp. Then the loop integral is divergent when k is parallel to xp. Using dimensional 
regularization, this collinear divergence appears as 1/e pole. These divergences may be 
factorized in the perturbative sense 

Mr = cr Pi , (51) 

where Cf^ is the finite coefficient function and Pi contains only the 1/e poles. On the other 
hand, the infrared-finite quantities F^b contain ultraviolet divergences which also show up 
as 1/e poles. When the infrared poles in Pj cancel all the ultraviolet poles in Fj^, T^^" is said 
to be factorizable. The product PiV-iB defines the renormalizcd parton correlation functions 
Fj. The final factorization formula for the Compton amplitude is then 

j.M-^^cT^r, , (52) 

i 

where Cf ^ is a well-defined perturbation series in ag and Fj is a finite nonperturbative 
distribution. 

The Fcynman diagrams in the hard part are simpler to evaluate after the collinear ex- 
pansion because all the external momenta become collinear (p and q). The only transverse 
momentum components come from the loop momenta. Therefore, it is convenient to intro- 
duce the light-cone coordinates 

= ^ (ifc° ± k"") ; k± = {k\e) . (53) 
v2 

The momentum-integral now becomes d'^k = dk^dk^d^k±. The dk~ integral can be done 
using the residue theorem by considering the complex plane of k~. After this is done, the 
transverse k± integrals become straightforward. The remaining A;+ integral is, in general, a 
nontrivial multidimensional integral containing both the evolution kernels Pi and the finite 
coefficient functions Cf Although this integral can be quite complicated, its dependence on 
e has been made explicit by performing the transverse integrals. This allows us to separate 
the process of verifying factorizability from the calculation of the coefficient functions. 



E. Summary 

We are now in a position to write down a set of rules for calculating the Compton 
amplitudes at an arbitrary order, n, in 1/Q. First, one must calculate the hard scattering 
Compton amplitude for every combination of g 'good' field components, where g runs from 
2 to n + 2. Only diagrams in which q^^ flows through every loop are considered in this 
calculation. Initial and final state interactions such as those represented in Fig. 4 are taken 
into account via renormalization Z-factors and the inclusion of the 'bad' field components. 
The external momenta in these amplitudes are written as 
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(a) (b) (0) 

FIG. 4. Hard scattering diagrams which represent processes that occur either entirely before or 
after the hard scattering, (a) and (6) are taken care of by the inclusion of 'bad' field components 
along with 'good' ones, while the contribution from (c) is grouped into a renormalization Z-factor. 

kf = (kj ■ n)p'' + \^{kj ■ p)n^' + Xkf , (54) 

and the amplitudes are expanded about A = 0. The coefficient of A"^^~^ is convoluted with 
the correlation function, and factors of kj- and kj are eliminated in favor of partial derivatives 
acting on fields. Second, one considers processes which involve one 'bad' field component 
and anywhere from one to n 'good' field components. The above process is repeated, but 
here the relevant coefficient comes at order \^+'^~9-'i _ Next, amplitudes with two 'bad' 
field components must be considered. This process continues until we have exhausted all 
possible ways to obtain a suppression of (Aqcd/Q)"- In general, one must consider a coUinear 
expansion to order _\"+2-ff-2b g^jj amplitudes containing g external 'good' field components 
and h external 'bad' field components for which n + 2 — (7 — 26 is positive or zero. 

As exemplefied by Eqs. (0,^), we could choose to eliminate the 'bad' components in 
the beginning of the calculation. In this scheme, some diagrams representing processes that 
occur entirely either before or after the hard scattering, such as Figs. 4(a) and (6) (but not 
(c)), are reinstated. Rules for the onshell couplings and propagators that appear in these 
special diagrams are derived directly from the equations of motion. This procedure is used 
in Refs. . Throughout the rest of this paper, however, we opt for the former method. 



The dependence of on the nucleon mass M through ^ (Eqs. (|T0|,|T3|)) prompts us to 
replace M — > AM in C, before each expansion. This produces additional suppression of 
certain correlations, but does not significantly alter the above analysis because its effects 
can always be grouped into the physical vector q^. If quark mass effects are desired, one has 
only to replace ip^ with — irtq in Eq. (|3^) , consider irtq — > Xirtq in the amplitudes, and 
expand as before. 

We note here that these expansions may be asymptotic in nature and cannot be used 
to predict experiment to arbitrary precision because of the problems presented by infrared 
renormalons fEOll- 



IV. ONE-LOOP CORRECTIONS TO THE TRANSVERSE COMPTON 

AMPLITUDE St 

Armed with the formalism presented in the last section, we are ready to study the main 
subject of the paper: one-loop corrections to the qt scaling function. As we have explained 
before, it is the Compton amplitude St that is a natural object to calculate in Feynman- 
Dyson perturbation theory; we focus on it exclusively throughout this section. Part of the 
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results here have already appeared in our previous publications [^JT^ • More details are given 
here for the reader's convenience. 



A. One-Loop Corrections in the Non-Singlet Sector 

The non-singlet sector requires at least two quark fields, so the correlation functions we 
need to consider are exactly the same as those in Section 3C. Before we begin, however, 
some comments are in order. 

As with most quantum loop-corrections, our calculation contains divergences. As ex- 
plained above Eq. (^l]), these divergences are of infrared origin after usual renormalization 
and can be combined with the bare distributions to generate a finite, scale-dependent renor- 
malized distribution. Working in dimensional regularization, one can isolate these diver- 
gences and obtain a relation between the renormalized and bare distributions. This relation 
allows us to determine the evolution of Kq{x, y, Q^) and check the validity of factoriza- 
tion for our process. Of course, in order to use dimensional regularization to calculate a 
spin-dependent quantity, one must specify a prescription for 75 and e^^'^^ . In this paper, we 
choose 't Hooft and Veltman's prescription . 



In principle, to obtain the full amplitude one must calculate all relevent Feynman dia- 
grams. However, loops associated exclusively with external legs can be grouped into renor- 
malization Z-factors for the external fields. The absence of a physical scale at the massless 
poles of the physical external propagators reduces these factors to unity. In addition, when 
both photon indices are left open, fundamental symmetries allow one to reduce the num- 
ber of diagrams that actually need to be calculated. The crossing symmetry, for example, 
allows one to cut independent diagrams in half by asserting symmetry under the replace- 
ment (yU z/, g ^ —q). In our case, the freedom to exchange /i and v is destroyed by our 
asymmetric treatment of these indices. However, we know that the part of the amplitude we 
search for is antisymmetric in /i and z/, so we can obtain the full result simply by calculating 
all diagrams with a certain ordering of photon insertions and subtracting (g — > — g) in the 
end. Although these diagrams are not all independent when photon indices are left open, 
they represent the minimal set of independent contributions to our calculation. 

By inspection, one can see that there are four diagrams associated with ((?,&) = (2,0), 
sixteen associated with (3,0), and eight with (1, 1). These diagrams can easily be grouped 
according to color factor. All of the diagrams associated with (2, 0) and (1, 1) are propor- 
tional to Cp = {N^ — l)/2N in SU(A^), while those representing (3, 0) also include the factors 
Ca = N and Cp — C^/2 = —1/2N. Routing the transverse momentum through the quark 
line and performing the collinear expansion 

-ff+-ffi^±-^+--- (55) 



before evaluation, we see that the four (2,0) graphs actually represent twelve contributions, 
each with an insertion of the operator ^_|_ — >■ i where the partial derivative acts 

on the external quark field. Each of these contributions corresponds to a similar graph in 
which the operator insertion is replaced by a gluon line. The only differences between these 
contributions are the change in momentum induced by a gluon absorption and the associated 
gauge group generator. Since the generator contributes only to the color factor, the (2, 0) 
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( 21 ) ( 22 ) ( 23 ) ( 24 ) 



FIG. 5. Diagrams contributing to the hard scattering coefficient M^^~^ relevant to St at one-loop 
order in the non-singlet sector. (1)-(12) are proportional to the casimir Cp, (13)-(20) to Cf — Ca/'^i 
and (21)-(24) to Ca- 

contribution is simply the CiT'-part of the (3,0) contribution evaluated at zero momentum 
transfer, i.e. x = y. This implies that the sum of (2, 0) and (3, 0) is obtained by replacing 
—gAaiiJ^n) with iDa^lJin) in the C^-part of (3,0). Hence, the Cir-part of our amplitude is 
already expressible in terms of KgB{x,y) alone. Since there are no partial derivative terms 
proportional to Ca, this part must be gauge-invariant by itself. Looking at Eq. (^Qf ), we see 
that the necessary requirement is that its coefficient vanish at least linearly a.s x ^ y. We 
will see below that this is in fact the case. 

Now that we have successfully reduced the number of diagrams we must calculate from 
fifty-six (actually, more like seventy-two) to the twenty-four shown in Fig. 5, we are ready 
to begin the actual calculation. Using the fact that /iip^ =^4'+ = ^ind that none of the 
external momenta have transverse components, one can reduce every numerator to simple 
Dirac structures times functions of k'j_, p ■ k, u, n ■ k, x, y, and Xb, where k^^ is the loop- 
momentum. Using techniques outlined in the Appendix, one can always cancel propagators 
until there are only three left. At that point, the integrals can easily be done by contour. 
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Useful integration formulae can be found in the Appendix; here, we simply give the result: 



.NS{1), 



2n 



Ns ( X y 



— , — ) -{xb ^ -xb) 
xb Xb 



KqB{x,y) ; (56) 



1-x 



+ 



'22/(1 -x) 2 
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log(l - x) - - log(l - y) 

X 1 — X J y 



Ay 
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+ 



+ 
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X 



x^ 



X 



x-y 



X . 



log(l — x) 
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x^ 



---^]log2(l-x) + -log2(l-2/) 
X 1 — X J y 



(57) 
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log(l - x) + 
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where = — (e^) is the nonsinglet part of the quark charge, (e^) = S^e^/n/ is the 
average squared quark charge for nj quarks. We have taken /i as a reference scale for our 
regularization procedure, and 7 = 0.57721 ... is the Euler constant. One can check explicitly 
that x-y divides the Ca part of this expression, as required by gauge invariance. 

Interpreting the divergent part of this result as the renormalization of our distribution, 
we find that 

«s(Q^) f dxdy 



X C, 



Xb x + xb 

y 



Xb — X \ 



^ ( M^B - X) ^Xb 
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(58) 



y \x-y xb{xb 
-Xb) \ KNs{x,y,Q'^) ■ 



Xb 



From this result, one can immediately see that the evolution of q^^ {x^Q"^) is not self- 
contained. This is one of the main reasons Eq. (H3) is misleading. The distribution 



Qt ix,Q ) makes sense by itself only for one value of Q . Once that value is changed, 
one requires input from the full function Kns{x, y, Q"^) to define it. 

At this point, we can take the imaginary part of Eq. (|56|) as xb approaches the real axis 
from below to obtain the next-to-leading expression for Qx^ixBiQ"^)- This is done in Ref. 
H]. For now, we turn to the singlet sector. 



B. One-Loop Corrections in the Singlet Sector 

At this order in ^^(Q^), gluons can interact with the photons through a virtual quark 
loop. This process generates a contribution to St from purely gluonic effects. Using our 
power-counting techniques, we see that once again we must consider three different cor- 
relations : = (2,0), (3,0), and (1,1). As before, the contributions to (2,0) can be 
obtained from those to (3,0). However, due to the bosonic symmetry of the gauge fields, 
the symmetries associated with this process are quite different. 

Let us consider first the contributions to (1, 1). These are represented by a quark loop 
with two photons, a transverse gluon, and a longitudinal gluon attached. Including the 
bosonic symmetries for both gluons and photons, a total of six diagrams contribute to 
this amplitude. Using the properties of fermionic traces, we can reduce this to the three 
contributions shown in Fig. 6. An explicit calculation shows that these contributions exactly 
cancel. Since this is the only place where, through Eq. (P^D, one can generate a purely 
singlet quark contribution, we see that the singlet quark contribution to S'>^'^ is identical to 
the nonsinglet contribution for each quark. 

The amplitude associated with (3, 0) has twenty-four separate diagramatic contributions 
once one has taken the three- fold bosonic symmetry into account. Leaving our gluon indices 
open, we can reduce this number to four. Once again, the properties of the fermion trace 
allow us to cut this number in half. Hence we are left with only two independent contri- 
butions to our amplitude. These contributions are represented by the diagrams in Fig. 7. 
To obtain the full amplitude, we simply double the charge-conjugation even part of these 
diagrams and add the gluon permutations. A symmetry factor of 1/6 must also be applied 
to associate the amplitude with our correlation function. 

After simplification, the contributions from Figs. 7 (a) and (6) have the form 



I 



{el) 



^Tf n ) i-9B)^ A,{x, y)gT9l'' + B,{x, y)gfgT + C,{x, y)g^''gf , (59) 
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FIG. 6. Diagrams relevant to St at one-loop order in the singlet sector involving one 'bad' 
^luon field component. The sum of the three contributions vanishes. 





FIG. 7. The two independent contributions to M^^f^^"' 



^^-'(3,0)b 



-Tpf 
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abc 



-9b) 



{el) 



(60) 



respectively, where /'^'"^ are the structure constants of the group, Tp = 1/2 is the generator 
normalization, and x = x/ xb-, y = u/xb represent the relative momentum fractions. Here, 
we have ignored the C-odd rf"'"^ contribution which cannot contribute due to Furry's theorem. 
The full amplitude becomes 
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FIG. 8. A diagrammatic representation of the consistency relations among loop functions re- 
sulting from our freedom to choose the routing of the transverse momentum. 



(X3gB)al3'y{x,y) 



+Bi{y - X, -x) + B2{y - x, -x) - Bi{x, x - y) - B2{x,x - y) (62) 
+Ci(x, y) + C2(x, y) - Ci{y - x, y) - C2{y - x,y) , 
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(63) 



where we have multiplied by a factor of two to take both possible directions of fermion flow 
into account and included a factor of six for gluon symmetry. Due to the bosonic symmetry 
of the gauge fields, all of the terms in Eq. ( |6l]) are identical; our result can be rewritten as 



2 u 



f dxdyff—,^) {TsgBmx,y) 

J \Xb Xr/ 



(64) 



This amplitude must combine with that associated with (2, 0) to make a fully gauge invariant 
contribution to s!j^\i',Q'^). 

The contributions to (2, 0) are calculated according to the prescription mentioned before. 
We send a small amount of transverse momentum into the quark loop through the incoming 
gluon and remove it through the outgoing. Internal propagators which carry this momentum 
are then broken successively by the operator ij^. Since we have a loop, we are given a choice 
of which way to route the momentum. Either choice must yield the same result, so we 
immediately have a set of relations among diagrams. As before, our new vertices differ only 
trivially from gluon insertions. Hence the relations between the (2, 0) diagrams translate 
immediately into consistency requirements on the functions Ai{x,y), Bi{x,y), and Ci{x,y). 
Six relations follow from the two diagramatic requirements shown in Fig. 8. They are 



Ai{0,x) + Bi{x,x) +C2(x,0) + Ci(x,0) = , 
Bi{0, x) + Ai{x, x) + A2{x, 0) + Ai{x, 0) = , 
Ci{0,x) + Ci{x,x) + B2{x,0) + Bi{x,0) = , 



(65) 
(66) 
(67) 
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xp.a.a xp,e,i 
(•) 



xp,a,a 



(b) 



FIG. 9. The two independent contributions to M'^^^^^ 



A^ix, x) - C2{-x, -x) + ^2(0, x) - C2(0, -x) = , 
B2ix, x) - B2{-x, -x) + ^2(0, x) - ^2(0, -x) =0 , 
C2{x, x) - A2{-x, -x) + ^2(0, x) - ^2(0, -x) = . 



(68) 
(69) 
(70) 



These relations provide a welcome check of our calculation, as all of the expressions involved 
are extremely complicated. A similar set of relations can be derived in the non-singlet 
sector by allowing the transverse momentum to flow through an internal gluon propagator. 
However, in this case the relations we would derive arc not as useful. 

Using exactly the same technique that allowed us to derive the consistency relations, we 
can obtain an expression for the full contribution to (2, 0). There are eight distinct diagrams, 
of which only two are independent. The diagrams in Fig. 9 (a) and (b) have the value 



(2,0)a 



m: 



(2,0)6 



{TPS'"') 



{el) 



-A,{x, x)gr9l' - B^{x, x)gfgT - C,{x, ~x)gT gf 
^2(0, x)g^J'gT + 52(0, ~x)gTgi' + ^^2(0, x)^r^f 



(71) 
(72) 



respectively. Adding the remaining diagrams and including a symmetry factor of 1/2, we 
find the full contribution to (2,0) : 



rplJ.+ (9) _ rp {eq) 

-'-{2,0) -^F 



V J 



dx 



Pi{-x)gT9± +P2{x)g±gT +Piix)g±gT (r2gB)7^a(2;) ; (73) 



fif} -fa 



IJ,-f af) 



Pi{x) = A-i{-x, -x) - Ci(x, x) + 6*2(0, x) - ^2(0, -x) , 
P2{x) = Bi{-x, -x) - Bi{x, x) - ^2(0, x) + ^2(0, -x) , 

'^"/"^ A^(0)i9;3^^(An) 



(r2gB)a/3j{x) = / ^ ^ 



PS 



PS 



(74) 
(75) 

(76) 



Since our amplitude is required to be symmetric under (a;^ — > —xb), the fact that P2{x) is 
an odd function of x implies immediately that it is zero. This allows us to write the full 
(2, 0) contribution as 



T'''+ (s) — OT 

-'(2,0) — ^-^F 



{el) 



dxp,{x){V2gBYl{x) . 



(77) 



To properly combine with the (3, 0) contribution, it would seem that we must have 
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Pi{x) 



(78) 



However, using the relations (|65| - |70D one can show that 



Pi[x) 



(79) 



At first glance, this is disastrous. It looks very much like we have made a mistake. On the 
other hand, we have been very careful to include all of the appropriate symmetry factors 
and counting. Going back through the analysis, one has no choice but to conclude that it is 
correct. The answer to our problem lies in a more subtle place. Looking at (T3gB)^aix,y), 
we see that it has the peculiar property 



dx (r 



3gB 



{x,y) = 



(80) 



Hence there is an ambiguity in our result. Any function g{y) which is independent of x 
can be added to f{x,y) to form an equally acceptable amphtude. The contributing part of 
f{x,y) is actually the modified function f{x,y) — f{0,y). The nature of this ambiguity is 
that it does not affect the (3, 0) contribution to our amplitude. However, it certainly will 
affect the (2, 0) part when one sets x = y. One can easily show that the modified function 
satisfies 



pi{x) = (/(x, x) - /(O, x)) , 

as required by gauge invariance. Defining the gauge invariant function 
1 



(81) 



-^xy 



PS 



(82) 



^2gB aix,y) + ^3gB aiX-iU) 
J AtX 477 \ 

= ie-^''^MSo.KgB{x,y) , 
and substituting the explicit expression for /(x, y), one arrives at the full gluon contribution: 



5«^(.,Q^) = ^ii^l^(eJ) 

ZTT 



dx dy 
xy{y - x) 



MM — , — 

Xb Xb 

- (xb - 



-Xb 



KgB{x,y) 



(83) 



M^{x,y) = — 
xy 




e/2 r 



[y — 2x)(l — x) + 'iy{y — x) — 6 — (y — x) ) log (1 — x 



X 



+ UAx - 3y)(l -y)+ Ax{y - x) - 8 - (y - x) j log (1 - y) 



+ 



(x + ?/)(2x-3?/) 



x-y 



(1 - (x - y)) - xy log (1 - (x - y)) 
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-3{y + 2 ) (1 _ x) log(l - x) 



X 



X 



+ \^3{Ax - 3y){l -y) + Ux{y - x) - 20 - {y - x) j log(l - y) 

Qx^ + xy-9y^ am n ^ aa 

H (1 - [x- y)) log (1 - (x - y)) 

x-y 

~2 (^-^ " '^^^^^ -^) + ^yiy - a;) - 6 ^ - x)^ log^ (1 - x) 
-i UAx - 3y){l -y)+ Ax{y - x) - 8 ^ - x) j log' (1 - y) 



1 ( {x + y){2x-3y) 

2 V ^-y 



(1 - (x - y)) - xy log' (1 - (x - y)) 



(84) 



which is identical to that in Ref. p6| and agrees with the results of Ref. after one has 
transformed into their basis. 

As before, the divergent part of this expression generates the evolution of the leading 
distribution. In this case, it contributes an extra term to the singlet quark distribution. 



Ks{x,y,Q^)^Y.K,{x,y^Q') ' 
that is not present in the nonsinglet sector : 

dx{ ^—] Q' :4g|(x,Q') = ^^^')^/^- f ^-^y 

X — xb X + xbJ dQ^ 



27r 



x'^y^(y — x) 



{y - 2x)(xb - x) + 3y{y - x) - 6 - x) log 1 

X / \ Xb 

(xbx \ f y 
(4x - 3y){xB -y) + 4:x{y - x) - 8 - x) log 1 
y J \ Xb 

( {x + y){2x - 3y) \^ ^~y 
+ I {xb - (x - y)) - x|/ log 1 



x-y 
■ {xb -Xb) 



Xb 



Kgix,y,Q^) + 



(85) 



(86) 



where the ellipses denote terms given by Eq. ( p8|) above. 

At this point, one can easily obtain an expression for QtIxb, <5') by taking the imaginary 



part of Eq. (|83|) as x^ approaches the positive x-axis from below. This is done explicitly in 
Ref. m. 



V. THE OPERATOR PRODUCT EXPANSION 

So far, our results have been expressed in terms of nonlocal parton correlation functions. 
While these distributions are easy to interpret, they do not allow one to exploit the full 
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symmetry of our theory. In particular, they do not lend themselves to a straightforward 
twist analysis. Furthermore, their evolution equations are unnecessarily complicated. 

To obtain an expression with explicit twist separation, we cast our result in the form of an 
operator product expansion of two electromagnetic currents. The idea behind this approach 
is that for large Q^, the electromagnetic currents in the definition of T'^'^ are separated 
approximately along the light-cone direction n^. Such operator products can be expressed 
as an infinite sum of local operators in the same way that a Taylor series can be used to 
represent a function. Since local operators can be separated into irreducible representations 
of the Lorentz group, they are easier to classify and their evolution is somewhat more 
constrained. 

One can obtain the operator product expansion by expanding our previous results about 
Xb = oo : 

ST{u,Q') = J2sn{Q'){ — ) . (87) 



n=0 



Xb/ 



Although this expansion takes place in the unphysical region xb > 1, its terms are related 
to moments of the physical structure function, qt^Xb^Q"^) through dispersion relations : 

dxx'''-gT{x,Q^) = -.s^{Q^) . (88) 
Jo 4 

This shows why the operator product expansion is experimentally relevant. 
At leading order in a^, it is quite trivial to expand our result. We obtain 

(g^) = 2j2el f dxx'-qT{x,Q') . (89) 

g "'-1 

From the definition of qrix, Q"^), one can show that 

rl 1 / _ 

J^^ dx x^" qrix, Q^) = [PS ipgiD^^iD^^ ■ ■ ■ iD^^^-i^-i^ip, 



(90) 



Substituting this result into Eq. (p8| ) gives us an expression for the (2n + l)-th moment of 
the physical structure function in terms of the matrix element of a local operator. 

The different components of the operator above do not constitute an irreducible repre- 
sentation of the Lorentz group. As such, they are not constrained to have simple properties 
as the renormalization scale is changed or higher-order radiative corrections are taken into 
account. A more useful basis is one which forms irreducible representations of the Lorentz 
group. Such a basis will be useful at arbitrary orders in perturbation theory, as the behavior 
of its elements under renormalization is restricted by the space-time symmetry. A standard 
analysis of the Lorentz group tells us that the components of the operators 

.Cm^' = ^,^D^'' ■ ■ ■ ^I^'^-T^"^^ VsV', , (91) 
.Cs = i'giD^'' ■ ■ ■ tD^'-Vh^^g , (92) 

where (• • ■) denotes symmetrization of indices and the removal of all traces and [■ ■ ■] indicates 
antisymmetrization of indices, form irreducible representations. Taking 



27 



PS 

PS 



PS) = 2<(Q2) MS'^^'P^^ ■ ■ ■ P^"+' 



(93) 
(94) 

we see that the matrix elements of every component of 6n,3 are kinematically suppressed by 
at least one power of Q in relation to the +• — h component of On+i,2- This suppression 
is characterized by the twist of the operator, defined as mass dimension minus spin. Since 
9n+i,2 has spin n + 1 and 9n,3 has spin n, the fully symmetric operators have twist-two and 
the operators of mixed symmetry have twist-three for all values of n. It is generally true 
that operators of twist t are kinematically suppressed by at least (Aqcd/Q)*~^- However, 
as exemplified by relations in Eqs. ( p3t |9^) , all of the relevant information can be extracted 
at 0(Aqcd/Q)*-'. 

In our case, St contains contributions from twist-two and -three. However, since the 
twist-two information can be extracted from the leading scaling function (7i(xb,Q^), only 
the twist-three contributions are new. To see this explicitly, we transform Eq. (PO) to the 



new basis and write [IS 



(95) 



An analysis of the leading contributions to T^'^ leads one to the moment relations 

dx X- g^'^ {x, Q') = Iy: el aUQ') , (96) 



which allows us to isolate the pure twist-three information [23]: 



n+1 



n 



2(n + 



(97) 



One of the main purposes of this paper is to present the one-loop corrections to Eq. 
In particular, it will be interesting to see whether or not the twist-three information is 
isolated by the same combination of qt and gi at higher orders. A priori, there doesn't seem 
to be any reason why it should be. This statement is really nothing more than a relation 
between gi and the twist-two part of gr- It seems to depend on the operator (^) relevant 
to gr- In fact, we will see that this particular combination of gr and gi does isolate pure 
twist-three effects at one-loop order. This fact asserts a kinematic rather than dynamic 
relationship between these two quantities, as will become clear below. 

In order to separate the twist-two and -three parts of our one-loop amplitudes, we must 
find an unambiguous way to tell the difference between these contributions. Since one 
cannot form a mixed-symmetry operator with only partial derivatives, d± cannot appear 
in the twist-three part of the expansion. This allows us to use the partial derivative terms 
in our distributions to identify the coefficients of the twist-two operators. Completing the 
operators using the definition in Eq. (0), we identify the remainder of the expansion as 
purely twist-three. 

Before we can systematically express the twist-three contribution as a sum of local oper- 
ators, we must identify an operator basis in which to express it. Unfortunately, the operator 
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6n,3 is not sufficient for our needs. One can immediately see this by looking at our expression 
(^) for the one-loop contribution. Our expansion will involve general moments of the func- 
tions Kq{x,y,Q'^) rather than the simple functions qT{x,Q'^). These moments are directly 
related to matrix elements of operators of the form 

ii^(id+)\Dl(id^y flj.i:, , (98) 

but these operators have no definite twist and are not useful for our purposes. A useful basis 
of twist-three operators in the nonsinglet sector was first identified in |20]. The operators 



1 



QCrfll---IJ,n _|_ Q<Jfll---IJ,n _i_ TJCT fJ-1 ■ ■ ■ fJ-n TjCr fj,l ■ ■ ■ fin 

1^n,j "Tg ^n,n-j n,j <? ^n,n-j 



for j = 1 to — 1, satisfy 
dx dy 



xy 



where 



PS 



1 n,j 



PS) = 2wl AQ^)MS^''P^'''^ ■■■P' 



(99) 
(100) 
(101) 



(102) 



(103) 



Note the appearance of the field strength tensor in our operator basis. Since there can be 
no partial derivative contributions to the twist-three part of our expansion, this is required 
by gauge invariance. 

Expanding S^^^^^ in Eq. (|56| ) about = oo and using Eq. ( |102| ), we find 



n even=2 



Xb 



n 



2NS 



^3NS 



CFUTl{n)-AS2{n) + 



2 2 

+ 



n+l n+2 



n-l 

-E 

Si(n) + 



^3NS ,„q 



(n - j) 



2 



n + 1 



n 



^2NS 



n+2 n+1 
Ca 1 



(104) 
(105) 



2(n + 11 



2 n + 2 



(2 + 5i(n + l)) 



(106) 



+ (-1)^ 



+ 



n 



(2 + 5i(n + l)) 



'n — \ 



n -|- 1 \ j — 1 



g^A;-l\2 + 5i(A;) 



k=3 



k 



where we have defined 



29 



Tf{n) = 



1 

i=i 

Sj{n) 



(107) 
(108) 



1=1 



In the above form, it is easy to see how much the result deviates from the naive ex- 
pectation that only the twist-three part of qT{x,Q^) will contribute to 5*^. The 'simple' 
twist-three operator 6^ 3 takes the form 



1 

An 



n — 1 

E< 



n 



(109) 



in the complete basis |2^. In deriving this result, the equation of motion has been repeatedly 
used. Also, like many of the operator relations in this paper, this equation is valid only for 
forward matrix elements. Looking at Eq. (|106|), we see that only the first term respects the 



naive expectation. Taking the large- A^c limit, which simplifies the evolution of qrix, Q ) as 
we will see below, gets rid of most of the complicated j-dependence of our result, but fails 
to remove the second term. All of the w^j enter nontrivially in the NLO result for S^^. 
This fact will complicate precision analysis of grix, Q^) data immensely. 

Since the operators Wnj for different values of j all have the same mass-dimension and 
transform in the same way under Lorentz rotations, they may mix with each other under 
renormalization. This is the origin of all of our difficulties with the evolution of qT{x,Q'^). 
Since we have no more symmetries to constrain the form of our result, a different combination 
of the Wnj will appear at each new order in a^. 

Expanding Eq. (|58|) about xb = 00, imposing order-by-order equality, and transforming 
to the standard basis, one arrives at evolution equations for the scalar matrix elements of 
our local operators : 



27r 



3 1 

2 (n + l)(n + 2) ' 



1 



3 

2 ^ (V^+l)(n + 2) 



+ -1 




-2Si(n + l) 



;iii) 



(^-jXf(Q^) • 



From these equations, it is quite clear that diagonal evolution of Ns^n,3 is broken only by 
terms of order in the large- A^c limit. This fact was first recognized by Ali, Braun and 



Hiller |2T| in 1991. Its anomalous dimension in this limit is reproduced by this equation. The 
evolution of the w^^iQ"^) was first obtained by Bukhostov, Kuraev, and Lipatov (BKL) pO| , 
and our result agrees with theirs. Therefore, we conclude that the amplitude is factorizable 
at this order. In addition, we reproduce the well-known anomalous dimensions of the twist- 
two operators Ns(^r, 



-1,2 
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In the singlet sector, things become somewhat more comphcated. Here, we must intro- 
duce an extra factor oil/xy to obtain gauge-invariance. To generate a local gauge-invariant 
operator product expansion, it would seem that this factor must be canceled by the kernel 
M^{x,y). As argued above, one must define M^{x,y) in such a way that M^(0,y) — 0. 
Hence, the only nontrivial requirement seems to be M^(x, 0) = 0. Checking this exphcitly, 
we see that it is not the case. Prom this result, one might conclude that a local gauge- 
invariant operator product expansion does not exist for St in the singlet channel. However, 
a more careful analysis brings a quite different result. 

The first step in a systematic approach to this problem lies in identifying and separating 
out the twist- two part of our amplitude. Since we know that partial derivative terms belong 
exclusively to this part of the amplitude, we can use them as a guide once more. In analogy 
with the nonsinglet sector, we write the twist-two operators 



(112) 



where the dual field strength tensor is defined by F'^" = ^e^'^^^Fa/s- Taking the (i -|- • • • -|-) 
component and simplifying, we obtain 



n+1 



[n + l)iA^{id+)''id''A'' + iA^{id+)''-^DaF''+ 



n—l 



+ P^''A^^{id+Y-^gAlid^Al + /"^"A^,(i(9+)™^A^(i(9+)"-'"A^ 

m=0 



(113) 



Here, the indices {i,j,k) represent transverse dimensions. Note that a sign has been taken 
into account for transverse partial derivatives and gauge fields with contravariant indices. 
There are several important features of this operator. One is the appearence of the operator 
DaF°'~^ . Through the gluon equations of motion, this term transmutes into a quark singlet 
contribution. It will generate a difference between the behavior of the singlet and nonsinglet 
twist-three quark operators under renormalization and their contributions to the operator 
product expansion. At present, its most striking feature is that it contains terms which can 
be generated by A'P{x, 0). This suggests that separating the twist-two contribution, whose 
gauge invariance and locality are not in question, from the amplitude before worrying about 
the issues of gauge invariance may lead us to the correct understanding of our result. 

In order to perform a straightforward separation, it is advantageous for us to write the 
three-gluon part of our amplitude in such a way that the bosonic symmetry is evident. As 
we have seen above, forgetting this symmetry can easily lead to incorrect conclusions. To 
this end, we define the local operator 



and re-write our twist-two operator in terms of it : 



(114) 



/)«+■•■+ 



^^e'^ |(n + l)iAHD^{id+)''Al - iA^{id+Y-^D^F'^+ 



n-l 
k=l 



n — l 



+ 



n 



+ 1> 



(115) 
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This equation is valid modulo total derivatives. Note that bosonic symmetry implies O', 
—On^n~j- Our full amplitude takes the form 

^+(1) _ asUfTF {el) - 



X {A 



11 even=2 



n+2 n+ 1 
2 



AHDtiid^TA", 



P5 



n-l 



— [^i(n) + 2]y, 
n + l)(n + 2) ^ ' ^fr{ + 



PS 



C ■ 



PS 



(116) 



where the divergent part has been ignored. It can be restored via the substitution 

2 /Q2g7\-^/2 



SAn) 



+ SUn) 



(117) 



In this form, it is easy to separate the twist-two and -three parts of the amplitude. We 
have only to identify an appropriate basis of twist-three gluonic operators to write the final 
result. An appropriate basis takes 
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for j = 2 to n — 1. The analogue of qRn,j is not independent in this case since 

F''+{id+YF^+{id+yF\ = F''+{id+yF\{id+yF^+ + F^+{id+YF''+{id^yF\ . (120) 

Note that our basis differs from that found in Ref. [Q. Through repeated partial integration 
and binomial coefficient sums, one can show that the relation 



n-l 



n— 1 -I n— 1 



A+-+ 



;i2i) 



j=2 ^ k=j 



is valid for forward matrix elements. It is this expression that establishes the gauge invariance 
of our result. Defining scalar matrix elements in exact analogy to the nonsinglet sector, we 
find the following result for the singlet part of s!^^'^ : 



2tt 



I) E 

n even=2 



2n 



(n + l)2(n + 2) 



-2<(q^)+e(":i^ 



1 

+7^ 



1 



n — 1 



2 (n+l)(n + 2) A.^nj^n,,l'^ 



neven=2 ^"^B 



in 



n-l 



n 



(122) 
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n — 1 



+ 1 



+2(-l) 



n + l 
n — 1 



J 



^ ^^31 — ("^iH + i) 



Here, the scalar matrix elements of singlet quark operators are defined via 

af(Q') = E<W') ' 



(123) 

(124) 
(125) 



This result was presented in a different operator basis in for n = 2, 4, 6. 

The evolution in the singlet sector can also be deduced by our results. One obtains 
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2 ^ {n + l){n + 2) 



2nfTp 
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-2Siin + l) 



(126) 
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2 (n + l)(n + 2) 
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n+l 



+ 2( 



^ ./^- 1^ (^ - 1)' + 2jn 
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(127) 

n-j)<.(g2) 



Here, also, one reproduces the well-known twist-two result p[ . Use of the equation of motion 
in the gluonic twist-two operator has generated a term which destroys the autonomous 
evolution of s(^n,3 in the large- iVc limit. This behavior can be predicted simply by looking 
at the diagrams which contribute to evolution in this sector, c.f. [^]. These results were 
first obtained by BKL |2D]. Translating from their basis to ours, we find that there is some 
disagreement between our results and theirs. 

We are now in a position to present the full NLO corrections to Eq. (P7\). Using the 
well-known result [^3H13 



x"^i(x,Q^ 



4tt 



^'''i' [n + l){n + 2) 



27r 



{S,{n) + l)ai{Q') 



(128) 
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one finds 



n+l 



gi{x,Q' 



E 



_ a.(Q^) ^ + 1 37VS 
27r n-j 



8(n+l) ^ 
(^-J)<,(Q') 
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2n 



vr (n + l)(n + 2) 



+ 1) E _ i)<AQ') (129) 



+ 



^^^^ E <Am I + ^ (aqcd/q, «^(q^ 



As promised, the relation between gi and the twist-two part of qt is unaffected at this order. 

We can gain insight into this phenomenon by stripping away the external states implicit 
in Eqns.( p!0^jl22| ) and writing the operator product expansion in its full glory : 

^ neven=2 W ) 
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as{Q^)nfTF 
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(130) 



a,(Q^)n/Tj. 
47r 



1 



ji-i 



(n+ l)(n + 2) 



i=2 



In this form, it is easy to see that the relationship between g\{x, Q"^) and the twist-two part 
of grix-, Q"^) is kinematical in nature. The operators involved are identical, so one expects 
their dynamical coefficients in this expansion to be identical. Stated another way, the twist- 
two contribution to this expansion is independent of the kinematics associated with the 
matrix elements we choose to take. From this point of view, one expects this behavior at 
all orders of perturbation theory. Inverting the Mellin transform, one obtains the functional 
relationship 



gi{x, Q') = g'^''{x, Q')x6{x - 1) - 9t-\x, Q') 



131) 



which must be satisfied if the operator product expansion is truly well-defined. 
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VI. FUNCTIONAL TWIST SEPARATION 



Although we have managed twist separation for the individual moments of the structure 
function, it is natural to ask if this separation can be done in term of the correlation functions 
themselves. This procedure would allow us to express fi'T*~^(a;B, Q^) in terms of twist-two 
distributions, for example. The most straightforward way of achieving this involves a simple 
resummation of the separate contributions to Eqs. ( |104| , |122| ). However, a quick glance at 
their form is more than enough to motivate us to find another way. 

The twist associated with a certain functional form is only well-defined after expansion 
since, as mentioned above, nonlocal operators do not form irreducible representations of the 
Lorentz group. However, it is quite easy to construct functions whose expansions display 
the symmetries of a given twist. For example, the form of the non-singlet twist-two operator 
(pT|) implies that for any function f{x) analytic in a region about the origin, the function 



T2ix,y) = f{x) + 



xf{x) _ yfiy) 
X — y X — y 



(132) 



will lead exclusively to twist-two local operators upon convolution with Kg{x,y,Q'^). Simi- 
larly, it is obvious from our operator basis (|101|) that for any function f{x, y) that is regular 
at X = the function 



Tz{x,y) = iy-x)f{x,y) 



(133) 



will lead exclusively to twist-three local operators upon convolution with Kg{x, y, Q"^). This 
fact is especially useful since it allows us to isolate the part of our amplitude which does not 
vanish as x ^ y into a function only of x, with corrections that contribute directly to the 
twist-three kernel : 



f{x,y) = limf{x,y) + 



f{x,y) - limf{x,y) 



(134) 



Now, we need only construct a way to separate the twist-two and -three contributions to a 
function only of x. 

In the non-singlet sector, this separation is quite straightforward. Since convoluting pure 
powers of x with Kg{x, y, Q"^) generates moments of qrix, Q^), one finds that the expression 



X 



n-l 



n+1 



n-l 



X 



n-l 



n—i—l 



j=0 



2n 1 
n + 1 2n 



n—l 



nx 



n-l 



- y: x'y''-'-' 



i=0 



(135) 



separates the twist-two and -three components of ^. 
functional relationship quite easily : 



This can be translated into a 



+ 



g{x) + 
d 



xg{x) 



ygiv) 



X — y X — y 

xg{x) yg{y) 



X — y X — y 



1 

^ ^ Jq dx' x' f{x') 



(136) 
(137) 
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Hence for any kernel M^^{x,y) in the non-singlet sector, the twist-two and -three parts are 
given by 



x-y 



x-y 



d A^<?/ X xm^^ix) ym^^(y) 

— xm^^(x) — + — + 

ax X — y X — y 



M'^^{x,y) - lim M'^^x,y) 
y 



respectively, where 



X^ Jo y^x' 



(138) 
(139) 

(140) 



While conceptually simple, this functional separation is practically quite complicated. 
As an example, we examine the leading order. Here, 



dx' 



x^ Jo x' — 1 



^log(l-a;) , 



:i4ii 



which leads to 



dxdy 



x^ 



log 1 - 
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Xb 
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dxdy 
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y{x- 
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1 



log 1 - 
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Xb 

y_ 

Xb 



(142) 



Kg{x,y,Q'^ 



x{x — y) 



- — log ( 1 - 



X 

Xb 



X 

Xb 



(143) 



2 , A y 
— log 1 

y[x -y) \ Xb 



Kg{x,y,Q^ 



Obviously, one can separate the kernel (0) in a similar way. However, the result of such a 
manipulation is extremely complicated and not very illuminating; it will not be presented 
here. 

The advantage of performing this separation is that we can explicitly replace Kq{x, y, Q'^) 
with the simple distribution 



Ag(a;,Q2) 
m our expression for brp 



d\ 



An 

AfS{0)tw-2/ ^2 



e^^^' [PS 
i^,Q^). Since 



ip JO) ^-f5ipg{\n] 



PS 
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n-l 



+ J2x'y--'- 
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n + 1 
' AM 



PS 



\<{Q') = \j dxx^Aq{x,Q') , (145) 
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one can write 



- > e: 



r dx , 




/ — log 


{'--) 


/ X 


\ Xb/ 



)Aqix,Q') . 



(146) 



More generally, for any kernel M'^^[x,y) in the non-singlet sector, the twist-two part can 
be written 



- j dxxm^^{x)Aq{x) 



(147) 



where m^^{x) is given above. 

In the singlet sector, the separation is spoiled by the presence of DaF""^ in Eq. ( |113| ). This 
makes it impossible for us to express the twist-two gluon operator in terms of Kg{x, y, Q^). 
One can circumnavigate this problem by introducing the singlet quark distribution. Exam- 
ining the twist-two gluon operator (|113|) , we see that 



dx dy 
xy 



n-l 

xy'^-' + Y.y'^'"' 

i=0 



Kg{x,y,Q'^ 



+ J dxdyiy-x)-''Ksix,y,Q') = ^a-^iQ') (148) 
for n even. In analogy with the non-singlet sector, we use this to separate the amplitude 



into its twist-two, 



dxdy 1 f X y ^ , ^ 
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and -three. 



dx dy 
xy 
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parts. Here, we have defined 



y-x 

Xb 



X J y — X 
Ks{x,y,Q^) 
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Xb 



Kg{x,y,Q^) (151) 



m^ix) 



I dx' lim — ^ — MHx\ y) 
Jo y^x' y — x' 



and introduced 



Ag{x) 



d\ 



JXx /pg 



F"+(0)F+(An) 



PS 



(152) 



(153) 



X J An 

Once again, explicit expressions for these separate kernels could obviously be obtained from 
(p^). However, due to the form of that result, such manipulations are better left for numerical 
analysis. 
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VII. CONCLUSION 



We have presented a systematic way to obtain radiative corrections to the power- 
suppressed terms in inclusive deep-inelastic scattering. This method is quite general and 
can be applied to other factorizable processes such as deeply-virtual Compton scattering 
and Drell-Yan scattering with few changes. 

Using this method, we have calculated the next-to-leading order corrections to the Comp- 
ton amplitude Sxii'iQ'^)- This result allows us to derive an expression for the operator 
product expansion of two electromagnetic currents separated along the light-cone valid to 
twist-three and 0{as). Our results show that the Wandzura-Wilczek relation between gi 
and the twist-two part of Qt is valid to this order and the form of the expansion indicates its 
validity to all orders in perturbation theory. The large- A^^^ simplification of the evolution of 
qT^{x,Q^) does not extend to its coefficient function at next-to-leading order. The general 
distributions Kpfsix^y^Q"^) enter nontrivially even in this limit. 

It would be interesting to apply our method to higher twist processes. In particular, 
ignoring interactions in the nuclear wavefunction, certain nuclear twist-four matrix elements 
can be factorized into the product of two nucleonic twist-two distributions. In view of recent 
data from RHIC, corrections to these processes may allow us to clearly determine the changes 
in parton distributions induced by the nuclear medium. In addition, radiative corrections 
to the nucleon's unpolarized longitudinal structure function Fl{x, Q^) can be analyzed with 
this technique. 



ACKNOWLEDGMENTS 

The authors would like to thank A. Bclitsky, W. Lu and X. Song for collaboration at early 
stages of this work. The work is supported in part by the Director, Office of Science, Office 
of High Energy and Nuclear Physics, and by the Office of Basic Energy Sciences, Division of 
Nuclear Sciences of the U.S. Department of Energy under grant nos. DE-FG02-93ER-40762 
and DE-AC03-76SF-00098. 



APPENDIX A: 

This appendix is meant to provide some formulae we have found useful during this 
calculation. The Feynman integrals one comes across here are either of the form 



or 



/ 



{27rY k'^{k + yipY ■ ■ ■ {k + ynPY{k + x^p + qY ■ ■ - {k + x^p + 

dfik {2p ■ kY-\k\Y+'^-^ nLi(^ ■n + at) 

{2t:Y k^{k + yipf ■■■{k + ynPf{k + xip + qY---{k + x^P + qf 



(Al) 



(A2) 



where n-\-m-\-l < 5 and j < 3. However, the formulae we have derived are valid quite gener- 
ally; we require only that the parameters are nonnegative integers that lead to nonnegative 
exponents. 
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To perform these integrals, one may use the relations 

1 1 
2p ■ k = (k + xpY H (k + ypY , 



X — y 



y-x 

2 Z + X 



(k + xpV + ~ ' (fc + ypV , 
X — y y — X 



(A3) 
(A4) 



kl = i±^^(A; + xp + qf + i±^^(A; + yp + qf , (A5) 
X — y y — X 

where z = k ■ n, to systematically cancel propagators. With both of the above integrals, 
one can continue this process until there are only three propagators left. The remaining 
integral can be easily done by contour, leaving only one nontrivial one-dimensional integral. 
With a suitable change of variables, this last integral can always be performed with repeated 
applications of 



In order to express the value of these integrals, we define the functions 
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:9n+i{x, X -y, {ai}, x - xb) 
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Xb - (x- y)' 

where i runs from 1 to and P^{xj) is the symmetric product of order i of the n objects 
{xj}, i.e. 

(AlO) 
(All) 



Pg {Xi) = X1X2X3 + X2X3X4 + X3X4X1 + X4X1X2 . 

Note that gn is a symmetric function of its last n arguments and 

d 



9nix,y,{a^}) = 

dan+i 

With these definitions, flAll) becomes 



IGvr^ \ 47r/i 



9n+iix,y, {ai},a„+i) 



g2g7^^ '^/2 ^ m n fi+j+m~i{xk,yi>, {yi}, {Xk " Xb}, {as}) 



(A12) 



k=x i=\ rifc' {^k> - Xk) V^i< {ye' - ye) 

where the products are defined to exclude the singular point, {ye} represents i y^'s, {xk} 
represents m — 1 x^'s, and {ag} represents all j a's. The integral in ( [A^ ) is identical to this 
result, with hi+j+rn-i replacing /j+^+^.i. 

With the help of these formulae, the calculation presented in this paper can be done 
almost exclusively on the computer. 
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